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Abstract 

A theory for lifting equations of motion for charged particle dynamics, subject to given elec- 
tromagnetic like forces, up to a gauge- free system of coupled Hamiltonian Vlasov-Maxwell like 
equations is given. The theory provides very general expressions for the polarization and magneti- 
zation vector fields in terms of the particle dynamics description of matter. Thus, as is common in 
plasma physics, the particle dynamics replaces conventional constitutive relations for matter. Sev- 
eral examples are considered including the usual Vlasov-Maxwell theory, a guiding center kinetic 
theory, Vlasov-Maxwell theory with the inclusion of spin, and a Vlasov-Maxwell theory with the 



inclusion of Dirac's magnetic monopoles. All are shown to be Hamiltonian field theories and the 
1-^ . Jacobi identity is proven directly. 
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I. INTRODUCTION 



In conventional treatments of electricity and magnetism phenomenological susceptibilities 
are introduced to describe material media. Concomitant with the introduction of these 
susceptibilities is the idea that charge can be separated into bound and free components, 
current can be similarly decomposed, and based on these separations expressions for the 
polarization and magnetization of the medium are obtained. However, it is well-known to 
plasma physicists that such a simple characterization is not possible for plasmas, where 
particle orbits may transition from trapped to passing and, indeed, may exhibit complicated 
behavior that can only be described by the self-consistent treatment of the dynamics of 
both the particles and the fields. Because of these complications, tractable and reliable 
expressions for the polarization and magnetization are not so forthcoming, particularly when 
approximations are made and/or additional physics is added. 

The purpose of the present paper is to construct a general theory for the coupling of charge 
carrying particle dynamics, entities possibly with internal degrees of freedom described by 
a kinetic theory, coupled to electromagnetic-like field theories. A theory that is gauge- 
free and ultimately expressible without the introduction of vector and scalar potentials is 
constructed. Like Maxwell's equations, the Vlasov-Maxwell equations, and virtually every 
important system in physics, the theory will have Hamiltonian form. This Hamiltonian form 
will be noncanonical, following the program initiated in Refs. jj, [2 ] . 

The construction begins in Sec. [Ill with a set of ordinary differential equations that de- 
scribe the particle dynamics. This set of equations, which is the basic model of the matter 
under consideration, is assumed to have a very general Hamiltonian form, possibly with an 
unconventional phase space and with a Hamiltonian that depends on specified electromag- 
netic fields including the field variables E(x, t) and B(x, t), and possibly all their derivatives. 
The problem then is to lift this finite-dimensional dynamical system that describes the mat- 
ter to a gauge-free field theory with a kinetic component that is of Vlasov type coupled 
to an electromagnetic component of Maxwell type. The difficulty with this lifting program 
lies in the coupling of the two components of the field theory. It is shown in Sec. II I II that 
the construction given naturally results in a field theory that is also Hamiltonian. This 
assures that there is a consistency to the coupling. Because the Hamiltonian theory requires 
variational calculus, it is most convenient to discuss constitutive relations resulting from the 



matter system in this section as well. In Sec. IIVI several examples are presented, beginning 
with the usual Vlasov-Maxwell system, followed by a general guiding center kinetic theory, 
a theory that includes spin and, to show the generality of our construction, a theory with 
monopole charge where the Maxwell field is modified. Gaussian units are used for all exam- 
ples. Section |V] contains concluding remarks. In Appendix |X] of the paper there are several 
subsections with direct proofs of the Jacobi identity for Poisson brackets of the noncanonical 
Hamiltonian field theories. The first one describes an old calculation of the author that has 
not heretofore appeared in print, a calculation that contains several useful techniques. The 
other subsections contain analyses of the other brackets of the examples of Sec. IIVI 



II. A GENERAL ELECTROMAGNETIC KINETIC THEORY VIA LIFTING 

Consider a general dynamical system with an n-dimensional phase space with coordinates 
z = (z 1 , z 2 , . . . , z n ) and evolution determined by a Poisson bracket and Hamiltonian S as 
follows: 

z a = [z a ,£} = J ab —, a,b = l,2,...,n (1) 
where the Poisson bracket on phase space functions g and h is defined by 

[9,h} = ^J^ b (2) 
L 1 dz a dz b 

and repeated indices are to be summed. The only requirements placed on the cosymplectic 
tensor J is that it endow the Poisson bracket with the properties of antisymmetry, [g,h] = 
—[h, g), and the Jacobi identity, [g, [h, k]) + [h, [k, g\] + [k, [g, h]] = for all functions g, h, k 
(see e.g. j^l, ij]). In the context of geometrical mechanics this is referred to as a flow on a 
Poisson manifold, but this formalism and its language will not be used here. Rather, the 
physics of matter described by this finite-dimensional dynamical system, as embodied in the 
Hamiltonian function S, is emphasized. The description of the matter in the formalism of 
this paper is contained in this function £ and its associated Poisson bracket (j2J). 

Particle orbits in given E(x, t) and B(x, t) fields are usually described in terms of the 
electromagnetic potentials, </>(x, t) and A(x, t), where 

1 <9A 

E = -V0 — and B = V x A . (3) 

c at 



Following this usual procedure, the Hamiltonian £ of the general system of ([I]) will be 
restricted for the purposes of the present lift theory to have the following form: 



£ = JC (p - eA/c, w, E, B, VE, VB, ...) + &{>, (4) 

where K. is an arbitrary function of its arguments. This form was proposed in the context 
of the variational theory of Refs. (J 7|. Since the aim is to generalize usual charged particle 
dynamics, the parameter e that denotes charge is included and c is the speed of light as 
usual. This particular form assures electromagnetic gauge invariance. Here, the phase 
space has been split into two parts, z p = (x, p, w 1 , w 2 , . . . , w d ), where the coordinates (x, p) 
are the usual canonical six-dimensional phase space coordinates and the coordinates w = 
(w l , w 2 , . . . , w d ) describe additional degrees of freedom, such as might occur in the classical 
description of a molecule with rotational or vibrational degrees of freedom. These additional 
coordinates will be referred to as internal degrees of freedom. An important, but not the 
only, Poisson bracket is given by 

[f,g] = [f,g] P + [f,g]w 

dg df df g dg 010 a /k\ 

=:V/ -^- V ^ + ^ J » ^' «>P = M->*> (5) 

which is the sum of canonical and internal pieces, and is assumed to satisfy the Jacobi 
identity. Such a bracket and Hamiltonian generate dynamics of the form of (pQ). 

Several comments on the Hamiltonian form of (j3J) and (jSJ) are in order. Note that £ may 
depend explicitly on the fields, and all of their derivatives, and the same may be true for the 
tensor J of the Poisson bracket provided the Jacobi identity is satisfied for any choice of these 
fields. The Poisson bracket may also have explicit w dependence, but no direct coupling 
of the internal degrees of freedom to the fields has been made explicit. In general, such 
coupling would need to be consistent with the symmetries of interest for these variables. For 
example, if w were a spin variable, say s, then £ would depend on s • B. This case is treated 
as an example in Sec. IIV CI Also, observe that all of the dependence on the spatial variable 
x occurs through the fields. The omission of such explicit dependence in £ is appropriate 
for media where spatial homogeneity is broken only by the presence of the fields; however, 
the x dependence could be added for further generalization. 

Alternatively, a manifestly gauge invariant form is obtained in terms of the coordinates 
.,w d ), where (x, v) denotes the usual six- dimensional velocity phase 



space coordinates, with 



p := mv H — A , (6) 



and m denoting the mass of the charged particle. In terms of these variables (CQ) takes the 
form 

z v =[z v ,K] + -X d --E, (7) 

TIL 

where the bracket of ([5]) becomes the Littlejohn [8|] Poisson bracket in ((7j), 

[g,h] L = [g,h]v + [g,h] s , (8) 

where 

r , 1 1 / m <9o\ , r e „ / (9(7 <97A 

S,/i v = - Vs- — -V/i--f and ^,h B = — B- /x— , 9 
m \ ov ay / 77rc \av dv / 

is a (6 + d) x 3 matrix used to embed E into the force law, and 

K (v, w; E, B, VE, VB, . . . ) = K (p - eA/c, w, E, B, VE, VB, . . . ) • (10) 

Thus the electric field appears as an external force in addition to any dependence on it that 
may come through the function /C, and the electromagnetic potentials no longer appear in 
the dynamics. Note, in general E cannot be written as a gradient in order to combine it 
with the first term of (I71). The dynamics of ((Tj) with arbitrary Poisson bracket in terms of 
z v , possibly depending explicitly on z v , E, and B, can be taken as the starting point for 
lifting to a kinetic theory. 

Usual Lorentzian dynamics is given by /C = m|v| 2 /2: when K, is written in terms of p, 
the bracket of (jSJ) with £ yields the equations of motion for a particle of charge e and mass 
m subject to given electric and magnetic fields. Alternatively, the same equations are given 
from © with [ , ] L . 

Now, the finite degree-of-freedom system of (JTj) will be lifted. The first step is to lift 
the particle dynamics to a kinetic theory for determining a phase space density f(z,t) = 
f(x,v,w,t). This is easily achieved by the standard Liouville form 

^ + [/,/C] + -E-|f = 0, (11) 
at m av 

where the generalization to multiple species is straightforward. Clearly the characteristic 
equations of ( ITT]) correspond to the finite-dimensional matter model of ([7]). 



The second part of lifting is to describe the coupling to Maxwell's equations. This coupling 
is effected by introducing the energy functional 

K[E, B, /] := J dxdvdw /C / , (12) 

whence the following expressions for the charge and current densities are obtained: 

dvdwf-V- — (13) 
t, ^ f , , d/C £ d 5K „ 5K 

J(x,t) = ey^- / + -- + cVx-. (14) 

Inserting these expressions for the sources into the usual form of Maxwell's equations com- 
pletes the lift. 

From ( ITS"]) and ( I14p it is evident that the polarization, P, and magnetization, M, can be 
identified as 

5K 5K 

p(x ' t) = ~k and M(x ' t) = _ k' (15) 

which is consistent with the usual definitions of bound charge density, polarization current, 
and magnetization current, 

dP 

p 6 =-V-P, j p = -^T> and Jm = cVxM, (16) 

respectively. Although the manner of lifting embodied in pip . f fl3|) . and (|14p is straight- 
forward, because of the functional derivatives in ( [15]) the dependencies of P and M on the 
fields E and B may be very complicated and contain high order spatial derivatives. 

In Sec. IHIl it is shown that this manner of lifting results in a Hamiltonian field description 
of the coupled system. It should be emphasized that this construction does not require the 
explicit introduction of the vector and scalar potentials. However, from the Hamiltonian 
form using £ it is clear that it subsumes the description using 0(x, t) and A(x, t)). To see 
the explicit form, define the momentum phase space density by /(x, p,w,t) = /(x, v,w,t), 
which gives under the change v -h- p the governing kinetic equation 

f 0, (17) 
with [g, h] = [/, g] p + [f,g] w . The coupling to Maxwell's equation is essentially unchanged: 

/' & K 

dpdwf-V-— (18) 

J(x,t) = e j /^- / + -- + cVx-, (19) 
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as are the expressions for P and M. Using the chain rule expressions, 

9f 1 df „ - „ „ e df' A Bf df e df dA 

TT"^' V / = V/ TT'VA, and -± = -± ■ — , 20 

op m (TV mc ov ot at mc Ov ot 

it is not difficult to show that [g, h] transforms to [g, h] L and ( ITTj) transforms into ( TTTI) . 



III. HAMILTONIAN FORM AND CONSTITUTIVE RELATIONS 

Based on past experience, viz. Vlasov-Maxwell and guiding center kinetic theories, a 
natural choice for the Hamiltonian functional is the following: 

H[f, E, B] = K - |dxE ■ ^ + i- JdK (|E| 2 + |B| 2 ) 

= K + JdxE-P + ^- Jdx(\E\ 2 + |B| 2 ) (21) 



where P, as given by ffl5|) . is used as a shorthand in the second line, which one could rewrite 
in terms of D := E + 47rP. The Hamiltonian of (I2T]) is a generalization of the energy 



component of the energy-momentum tensor first derived by variational methods in [5|- 
It is straightforward to verify directly that (12 ip is conserved by the combined field theory, 
Maxwell's equations with the sources (IT5|) and (jHJ) coupled to the kinetic theory of f fTTj) . 

One might think that the |B| 2 term of (12 lj) should be replaced by B ■ H, where H = 
B — 47rM, but this is incorrect. All polarization and magnetization effects are modeled here 
by the terms involving /C, i.e., they are a consequence of the particle dynamics. Rather than 
relating E and B to D and H by constitutive relations, the particle dynamics, extended and 
other, describes the physics that is often approximated by simplistic constitutive relations. 
For example, the difference between |E| 2 and E • D arises from the K. term that contains the 
matter dynamical information. 

If only (/, E, B) are used as dynamical variables, there is a difficulty in obtaining a 
Poisson bracket description for the field theory. The problem is readily encountered when 
one attempts to include polarization effects, because the polarization current has a time 
derivative and Poisson bracket expressions such as {E, H} do not produce terms with time 
derivatives of the dynamical variables; i.e., in Hamiltonian theories all time derivatives are 
on the left hand side, so to speak. Consequently a term of the form dP/dt cannot appear. 
However, there is a way to circumvent this problem, a problem that does not occur in action 
principle formulations, such as those of js- 7| . 



Functional differentiation of ( 121]) gives 



5E 



\5E5EJ ' ' An 

5K ( 5 2 K \ 1 ^ B 
1 E -I 

47T 



5B ^5B8E ) ' E + 



(22) 
(23) 



where 5 2 K/5E5E and 5 2 K/5E5H are second functional derivative operators that satisfy 



[5 2 K\ ] 6 2 K , f 5 2 K \ ] 
[sESe) =6E6E and 



5 2 K 



(24) 



(See e.g. {3} 0, 9] for a review of functional differentiation.) Expressions f l22l and ( 1231) 
reveal how polarization and magnetization effects are embodied in /C. Since the functional 
derivatives above have, in a sense, 'dressed' E and B, existing Hamiltonian structures will not 
be adequate. It is clear that without some modification one cannot obtain the polarization 
current. 

If the theory were expressed in terms of D, the following bracket on functionals F\D, B, /] 
could give the correct temporal evolution of D, provided 5H/5H = H/Att: 

(25) 
(26) 
(27) 
(28) 

with [f,g] v and [f,g] B given by ®, d v g = dg/dv, and F f := 6F/6f^ D := 5F/5V, 
etc. The Born-Infeld term of ( f28i) is motivated by their original theory 10j that was also 



{F, G} = jd^dwdwf{[F f ,G f ] v + [Ff, Gf] w ) 

dxdvdw f [Ff, Gf\& 
Aire 



m 



J dxdvdw f (G D • d v F f - F D ■ d v G f ) 
+Anc /rfx (F D ■ V x G B - G D ■ V x F B ) 



written in terms of D and B (see also 



111]). Although this term can give something like 



dE)/dt = V x H, it remains to properly define the meaning of D and H. Thus, this bracket 
alone does not constitute a closed theory. Similarly the coupling term (127]) . a generalization 
of that introduced in {2], 3] that includes the internal variable w, is written here in terms 
of D, but the generalization of the Marsden-Weinstein term J12I (see also [k^) of ( 126|) . and 
the first term of ( 1251) . also a generalization of that given in J 3], are unchanged. The new 
internal term here of (1251) with [Ff,Gf] w does not depend on D and does not affect the 
Jacobi identity (cf. Appendix IA2j) . 



To close the theory requires a constitutive relation, something like D = e • E. Such 
relations are often appended to electromagnetic theory based on phenomenological material 
properties, but here they emerge as a consequence of the Vlasov-like dynamics and the 
definitions ( I15p . Using ffTBT) gives 

D = D[E, B; /] = E + 4ttP[E, B; /] (29) 

with both P and D linear in /, but not in E and B. In general these functionals can be 
nonlinear and even global in nature. It is only required that there be a unique inverse 

E = D- 1 [D,B;/] = E[D,B;/]. (30) 

Similarly, using ([15!) 

H = H[B, E; /] = B - 4ttM[B, E; /] , (31) 
which is also assumed to have an inverse, i.e. 

B = B[H, E; /] = H + 4ttM[H, E; /] . (32) 

For given /C, the expressions of (fT5j) can be quite complicated, particularly when deriva- 
tives of the fields are included. However, these expressions are local in time, i.e., E, B, D, 
and / are all evaluated at the same time. Because of the presence of / and the equation 
governing it, causal effects are included in a dynamical sense and do not need to be put in 
at the expense of breaking time-reversal symmetry. Also, there is no artificial separation 
of charge into bound and free components or current into magnetization or other. Rather, 
charges and currents are determined dynamically according to the Vlasov equation. A given 
charge may behave in any manner consistent with this dynamics. 

When limE-s-o D[E, B; /] = the first order term in an expansion in E gives D = e • 
E, where e is the dielectric permittivity operator. Similarly, B = jj, ■ H, where [i is the 
permeability operator. If one were to replace the Vlasov dynamics by trivial dynamics of 
linear response away from equilibrium, then one can recover the usual permittivity and 
permeability relations, including the usual causal (see e.g. [14]) form in space and time. 
But, this will not be done here. 

To sum up, the Hamiltonian of (1211) is given in terms of (E, B, /), the bracket of ( J25l) - (l28l) 
in terms of (D, B, /), and (129]) is a closure relation relating D to E. Thus, if H\D, B, /] = 
H[E, B, /] is defined by inserting the inverse of (T29l in for E, a closed theory is obtained. 
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For general K. this inversion cannot be done explicitly (although a series expansion may be 
possible); however, the chain rule can be used to relate functional derivatives of H to those 
of H and thereby obtain equations for the time derivatives of (D, B, /) which can then be 
shown to be equivalent to those of Sec. [TT1 Alternatively, the chain rule can be used to write 
the bracket of (j25p - (j28p in terms of the set of fundamental variables (E, B, /). 

To understand the chain rule, suppose two functionals are related by F[D,B,/] = 
F\E, B,/], where the right hand side is obtained by inserting (1291) in for D in the func- 
tional F on the left. Variation of F = F gives 

( 5F m 5F X t>\ L 6F ^ L ( 5F ~ 5F 5F 

• AD -A JR 4- It] 7 Sf = My I- 



dx \sd sd+ sb' sb ) + r jjV-rkw^w SB > + r w Sf • (33) 

while variation of (1301) gives 

6E 5E 5E 

SE = m- SD + M- SB + JJ Sf - (34) 

where 5E/5D etc. are the usual Frechet derivatives obtained by first variation. Inserting 
( I34jl into ( 1331) and comparing the coefficients of the independent variations 5D etc., gives 

5F _ f5E\ ] 5F 

5V-\5V)'6E> (35) 
5F 5F /<5E\ t 5F 



5B 5B ' \5Bj <5E' ^ 
SF = SF (SE\[SF 

Sf Sf \6f) 5E' 1 ; 

A more explicit expression for (5E/5D) can be obtained by varying (129 j) at fixed B and 
/, giving 

ro = (£-*sKi)-* B = : s- 4B < 38 > 

or 

SE 

5E = — ■ 5E> = e _1 - <5D (39) 
oD - 

where e is the nonlinear permittivity operator (not to be confused with e). Evidently, 



and Ed = (Ed)t or (s^ 1 )^ = It is important to note that although <5D = e ■ SE, 

D ^ e-E; the correct relation between D and E is given by the nonlinear expression of ( 1291 . 
Similarly, 

5E i S 2 K , _ 5E i S 2 K 

**--m =4 "*~ K mE and ~ ~&j ~ = ~ ■ («) 
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-i 



and the functional derivatives of (l35l) - (l3Tj) can now be calculated: 

-i . E 



Hy> = E D ■ H E — e ■ e • 

— — All 

E 

An 



V<5B5E J' + 4^ + V= 5B5E ) ' 



47T 



= /C . (42) 

Now, using the expressions of (J4"2"j) in the bracket of fl25l) - fl28|) gives 

= -4vrcV x # D = -cV x E (43) 
dt v ; 

<9D - 47re /" 

— — = Aire V x H B / d\dw fd v H f 

at m J 

Ane r 

= c V x H / dvdw fdylC (44) 

m J 

= -[/,/C]--E-^, (45) 
m dv 

where [, ] = [, )l + [, ] w as defined by (jSJ) and (jUj). Thus the bracket reproduces the 
Vlasov-like equation of (11 II) and Maxwell's equations with the polarization and magnetiza- 
tion currents, but remember / and [ , ] can be written in terms of p using ([6]) and thus the 
above is also equivalent to (1111) . 

In the above, D is a convenience, a shorthand for E + 47rP, with E being the fundamental 
variable. One could eliminate D from these equations e.g. by writing (jHI) as follows: 

<9D „ «9E „ «9B „ df 

lF = E -a +PB 'lir + p ^i- < 46) 

where Ed = e, as before, and Pb and P / are again operators obtained by variation of P, and 
then inserting the other two equations of motion for the time derivatives. This procedure 
will lead to a complicated set of equations in terms of the fundamental variables (E, B, /). 
Another way of obtaining these complicated equations is to obtain a bracket in terms of 
E,B, and / alone, by inserting the transformations for the functional derivatives of fl35|) . 
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(136]) . and (}3Tj) into the bracket of (|25|) - (|28|) . This yields the following complicated bracket: 
{F, G} = J dxdvdw f F f + E\-F E ,G f + EyG E 

Jd^dvdw / ((e{, ■ G E ) • d v (F f + E} • F E ) - (Ej, • F E ) ■ d v (G f + E+ ■ G E 
+4ttc Jdx ( (Ej, ■ F E ) ■ V x (G B + E+j- G E ) 

-(e^-Ge) -Vx (Fb + E^-Fe)), (47) 



Aire 
+ 

m 



where to complete the procedure the expressions of (1401) and (1411) are to be inserted. Thus 
the bracket of (T4TI) is quite complicated with many terms and operators. Nevertheless, by 
its construction it satisfies the Jacobi identity (modulo the V ■ B = obstruction discussed 
in Appendix I A II) . With Hamiltonian H [E, B, /] it is easily seen from (jl2l that this bracket 
produces the correct equations for df /dt and <9B/<9t, but it is less easy, yet possible, to see 
it produces the complications of ()46l) correctly. 

For some theories of interest, including usual linear response theory, K. has a simplified 
form, viz. 

/C(v, w, E, B) = h(v, w, B) + V(v, w, B) • E + ^ E ■ |(v, w, B) • E , (48) 

where = fc. For example, such a linear polarization theory is sufficient for some drift and 
gyrokinetic-like theories. With ( T48]) 



P = - Jdvdw^f = — J dvdw Vf — E • J iwdw }g f . (49) 

The first term of the second equality of (149!) represents a permanent dipole moment per 
unit volume, which can be dropped: because (|2T]) is a Legendre transform it will cancel out 
anyway. The second term of (l4"9l) defines the electric susceptibility, Xe- Thus P = Xe • E or 

D = |-E (50) 

with e = / + 47rx e . Note that unlike the e of (1381) . when (T48l) is assumed, e is independent 
of E but not of B. Explicitly in terms of indices 



€ ij 



6ij — An J dvdw fcy(v, w, B) / . (51) 



Assuming e 1 exists and has the form 



I" 1 = / + 4 ?rx = I - 4vr^ + (4vr|e) 2 - . . . , (52) 
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where Xe 2 stands for matrix multiplication, and so on down the line. Hence, E = e _1 • D. 
Although e is linear in /, this is not the case for e _1 . Finally, for the /C of fj48|) that is 
quadratic in E, the following simplified expressions are obtained: 

EL = I" 1 (53) 
t de- 1 

E B= D 'ir < 54 > 

Note, obtaining these formulas can be facilitated by using identities obtained by varying the 
expression e -1 ■ e = /. 



IV. EXAMPLES 

In this section four examples are given: that of Sec. IIV Al is the usual Vlasov-Maxwell 
theory, that of Sec. IIV Bl is a guiding center drift kinetic theory that includes nontrivial 
polarization and magnetization effects, that of Sec. IIV CI includes a physically perspicuous 
internal variable, and that of Sec. IIV Dl was chosen to show the generality of the lift theory 
by altering Maxwell's equations. 



A. Vlasov-Maxwell 



For Vlasov-Maxwell theory w is nonexistent and only z = (x, p) appears. Thus, /(x, p, t) 
and, with ]C = |p — eA/c| 2 /2m, Eq. ( TTTT) becomes 

|P - eA/c| 2 



df 

at 



2m 



+ e<j>J 



— (p - eA/c) ■ VA • ^ + eV0 • |^ - - (p - eA/c) ■ V/ . 
mc op op m 



(56) 



In terms of /(x, v,t), /C = m|v| 2 /2 and ( l56l) becomes 



W = _ v . v/ _±( E+ v xB v« 

ot m \ c / av 



(57) 



From the general Hamiltonian ( ETT]) . with /C = m|v| 2 /2, evidently the Vlasov-Maxwell 
Hamiltonian is 

H = jJ rfxdv |v| 2 / + i- |dx (|E| 2 + |B| 2 ) , (58) 
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and with [g, K\l the bracket of (!25!) - fT28l) becomes the Vlasov-Maxwell bracket 

f f Aire 
{F,G}= /rfxdv (f[F f ,G f ] v + f[F f ,Gf] B + — f(Gvd v F f -F E -d v G f , 

+4vrc Jdx (Fe ■ V x G B - G E ■ V x F B ) . (59) 

With this bracket and Hamiltonian, one obtains the usual Vlasov-Maxwell equation as 

| = {/ , ff} = - [/ , K]i -£ E .|, (60) 

which is equivalent to (1571) . Similarly, since for this example D = E, the usual expres- 
sion for the current is obtained from {E,H}, there being no polarization or magnetization 
contributions, and Faraday's law is given by dB/dt = {B, H}. 



The relativistic Vlasov-Maxwell theory similarly follows with the choice /C = cyfp — eA/c| 2 + m^c 2 
and the theory can be written in manifestly covariant form jt|, but this will not pursued 
further here. 



B. Guiding Center Drift Kinetic Theory 

A canonical Hamiltonian description for guiding center particle motion was obtained in 
Q, B] by applying Dirac constraint theory to Littlejohn's degenerate Lagrangian [jj], If]] 
(with a regularization suggested in [l^]) in order to effect a Legendre transformation. The 
canonical variables of the theory are (x, q^, p, mjand the particle Hamiltonian, Dirac's 



19), is 



primary Hamiltonian for this problem (see e.g. 

£(x, g 4 , p,p 4 , E, B, VB) = v s • (p - eA*/c) + V m + e0* (61) 

= /C(p-eA/c,g 4 ,P4,E,B,VB) + e0, (62) 

which is of the form of (jlj). Here 

A* = A + bmcwo7(g 4 /wo)/e, e<p* = e<p + //|B| + m (q\ + \w E \ 2 ) /2 , (63) 

v £ = c(ExB)/|B| 2 , b = B/|B|, (64) 

v 3 = g 4 B7(y^) + cE* x B/Bl , V 4 = eE* ■ B*/(mg'B$ , (65) 

1 BA* 

B* = V x A* and E* = -W<p* — . (66) 

Y c dt v ' 

The function 7(2) is an antisymmetric regularization function with z = q^/vo and Vq some 

constant velocity. The Littlejohn theory is recovered if 7(2:) = z, in which case g 4 = V\\. 
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In the regularized theory j(z) ~ z for small but for large 7 is bounded so that 
with t> >> Vthermah mvoc-y(oo) « e|B|/(b ■ V x b), which is accomplished, e.g., by 7 = 
tanh(z). 

This theory has an eight-dimensional phase space with the canonical bracket 

[9, h] = [g, h] p + d q4 g d P4 h - d q4 hd P4 g , (67) 



which with K of (162|) completes the theory. The appropriate bracket of the form of ([25]) - 
(|28p is obtained with (|67[) . and from K, the functional K can be constructed and thus the 
Hamiltonian of (|2"T|) . This bracket and Hamiltonian produces the equations of motion; thus 
this system is a Hamiltonian field theory. From K the polarization and magnetization are 



can 



)e obtained straightforwardly. Since expressions are complicated, the reader is referred 



to 16, 71 for details. 



C. Spin Vlasov-Maxwell 

The nonrelativistic spin Vlasov-Maxwell system is a kinetic theory gene ra 
Vlasov-Maxwell system that includes a semiclassical description of spin 
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ization of the 



21|. The Hamil- 



tonian description of this system [2lj] will be shown to fit within the gauge-free lifiting 
framework. The spin Vlasov-Maxwell electron distribution function, /(x, v, s, £), contains 
the internal spin variable s = (si, S2, S3), and satisfies 



df 
dt 



-V-V/ + 



^(E + !xB) + ^V(s.B) 

m V c / mfic 



where m and e > are the electron mass and charge, respectively, 2nh is Planck's constant, 
fi e = p/xb/2 is the electron magnetic moment in terms of /ie, the Bohr magneton, and the 
electron spin g-factor. Equation (1681) is coupled to the dynamical Maxwell's equations, 



dB 

~dt 



_ ^ dE 

■ cVxE - m 



cV x B - 4ttJ 



(69) 



through the current J = J/ + cV x M, which has "free" and spin magnetization parts: 



Jf := — e dv ds v/ 



and 



M :-- 



2/4 

h 



ds sf 



(70) 



Equations (!68|) and ( 169|) . with (!70|) . are to be viewed classically and consequently a full nine- 
dimensional phase space integration, d 9 z = dxdvcfe, is considered for /. Spin quantization 
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is obtained as an initial condition that constrains s to lie on a sphere (see 211]). Extension 
to multiple species is straightforward, so will not be included. 
For this system fC is chosen as follows: 

/C(v,s,B) = ^|v| 2 + ^s-B, (71) 

however, more general forms are possible. It is not difficult to check that the characteristic 
equations of fl68|) are of the form of ffTT]) with the Poisson bracket [g, h] = [g, h] v + [g, h] s , 
where 

\g,h] a = B-(d a gxd a h), (72) 



with d s := d/ds, and /C given by flTTT) . 

Thus the Hamiltonian functional (I2TT) becomes 



H[E,B,f] = fd*z (yM 2 + ^.B)/ + i- |rfx(|E| 2 + |B| 2 ) , 



(73) 



which can be shown to be conserved directly by using the equations of motion, and the 
bracket of ( |25l) -( |28l) adapted to the present example is 



{F,G} sVM = J <f>zf([F f ,G f ]„ + [F f ,G f ] B + [F f ,G f ) s (74) 

Attp \ 
+ (i*E • d v Gf — Ge • 9 v Ff) ) (75) 

+4vrc Jdx (F E • V x G B - G E • V x F B ) . (76) 

The last term of (1741) of { , } s vm accommodates the spin, an internal variable; it is not 
surprising that it has an inner bracket based on the so (3) algebra (e.g. [lij). The remaining 
terms of ([71]) . fl75]) . and (ITS"]) produce the usual terms of Vlasov-Maxwell theory: it is a 
straightforward exercise to show that Eqs. (168]) and (169]) are given as follows: 

d f (9B (9E 

= {f, H} sV M , -qT — {B, H} S VM , -qT = {E, ff}sVM • 

This is facilitated by the identity j d 9 zf[g,h] = — J d 9 zg[f,h], which works for each of 
[g, h] v , [g, h] B , and [g, h] s . It follows that the polarization P = and 5K/5~Q = — M as 
given by (TTD]) . 

In Appendix IA 21 a direct proof of the Jacobi identity for {/, H} s vm is given. 
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D. Monopole Vlasov-Maxwell 



This example differs from the previous ones in that Maxwell's equations are changed to 
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23| with monopole 



another field theory. In particular, Dirac's theory of electromagnetism 
charge will be treated in the lift framework. 

For Dirac's theory, (J7J) is replaced by the particle orbit equations 

x = v and v = — (E+ - x B) + — (B- - x E). (77) 

m c m c 

where g and e are magnetic and electric charges, respectively. The appropriate particle 
Poisson bracket for insertion into (p5l) - ([28l) is 

[g, h} m = [g, h] v + [g, h] B + [g, h] B (78) 

where 

[f,g]n '■= ^-E • (dyf x d v g) . (79) 

mrc 

The particle Hamiltonian is given by K. = m|v| 2 /2, just as for Vlasov-Maxwell theory. 

In lifting to a kinetic theory there are various kinds of multi-species dynamics, with and 
without magnetic charge, that could be considered. Here the case of a single species of 
identical particles that carry both magnetic and electric charges will be developed, along 



24). Inserting ([TED into dSSJ) ([28]) and 



the lines of the quantum fluid theory considered in 
adding a new coupling term to account for the B that acts as an external force, gives the 
bracket 



{F,G} mVM = J dxdv{f[F f ,G f } m (80) 

+— / (G E • d v F f -F E -d v G f ) + ^-f (G B ■ d v F f - F B ■ d v G f ) >81) 
m m / 

+4vrc J chc (-Fe ■ V x G B — Ge • V x F B ) ■ (82) 

With Hamiltonian of ( 1581) this bracket yields 

| = _ v . V/ _g.^ (E+ I x B) + 9(B-IxE)) (83) 
at av \m c m c / 

dE 

— = cVxB-4vrJ e (84) 
dt v 1 

dB , N 

— = -cV x E-4vrJ m (85) 
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where 

J e = e J fvdx and J m = g J fvdx. (86) 

Thus monopole Vlasov-Maxwell is a Hamiltonian field theory. 

One reason for investigating Dirac's model in the present context is to see if the V • B = 
obstruction to Jacobi discussed in Appendix [A] can be removed. In Appendix I A 31 the Jacobi 
identity for {F, G} m vM is proved directly, and there it is discovered that the solenoidal 
character of B is replaced by V • (eB — gE) = 0. Thus, the space of functionals must still 
be restricted to such fields, as discussed in Append ix IA1 However, Dirac constraint theory 
can reduce this to a boundary condition at infinity [25 1. 

Another reason for investigating monopole theories is for their utility in developing numer- 
ical algorithms. For example, the Gudunov numerical method for magnetohydrodynamics 
{^(J (see also j^) exploits a form that allows for V ■ B ^ that was subsequently shown 



28 



30|). However, 



to be Hamiltonian without the V • B = constraint in [l|, |3j (see also 
since the monopole theory of this section requires a specific linear combination of E and B 
to be divergence free, adaptations of these methods in not so straightforward. In any event, 
the reader can rest assured if mononpoles are discovered there still will exist Hamiltonian 
guiding center and gyrokinetic kinetic theories, obtained with suitable choices for /C with 
associated generalized polarization and magnetization vectors. 

V. CONCLUSIONS 

The main accomplishment of this work is to describe how a matter model of dynamics 
can be lifted to a Hamiltonian coupled Vlasov-Maxwell system. En route to this Hamilto- 
nian theory, the general constitutive relations of fl30|) and fl3Tj) or, equivalently, the nonlinear 
permittivity and permeability operators, as determined by ( jl5]) . were obtained. These consti- 
tutive functionals are very general: as discussed in the paper, K, may contain all derivatives 
of the fields and may even be global in nature and contain integral operators. From the 
general constitutive functionals, it was shown how to obtain the usual linear relations. A 
noncanonical Poisson bracket, another step in the program started in jl-3|, was obtained 
for this general class of theories. 

Four examples were given, including the general class of guiding center kinetic theories of 
Sec. IIVB1 This latter example, like all systems that are in the class of variational theories of 
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6), easily was shown to possess the Hamiltonian structure. Thus, the theory of this paper 
determines the path to follow for obtaining the Hamiltonian formulation of a consistent 



gyrokinetic theory by making use of the results |3ll . [32]. 

Various generalizations are possible. Namely, the extension to many species of different 
dynamics, relativistic theory other than the Vlasov-Maxwell example that was described, 
versions where particle matter models have more general finite-dimensional Poisson brackets, 
are all straightforward. Also, extending the matter dynamics by coupling to other general 
gauge-free field theories is possible. 

Another application of the techniques of this paper will be used in a subsequent work [33] 
that will treat Hamiltonian perturbation theory in the field theory context. There it will 
be shown how an exact transformation of the particle (characteristic) equations of Vlasov- 
Maxwell equations can be lifted to the kinetic and Maxwell equations, and how this can be 
used in perturbation theory for infinite-dimensional noncanonical Hamiltonian systems. 

In physics there are two ways of constructing new theories. The usual way is to con- 
struct an action principle by postulating a Lagrangian density with the desired observables 
and symmetry group properties. Alternatively one can postulate an energy functional and 
Poisson bracket, which is essentially the approach of the present paper. With this lat- 
ter approach, one must prove directly the Jacobi identity {{F,G}, H} + {{G, H},F} + 
{{H, F}, G} = 0, for all functionals F, G, and H. Techniques for doing this are not gener- 
ally known, and this provides one reason for Appendix [A] 



Appendix A: Direct proof of Jacobi identities 



One term of the original Vlasov-Maxwell Hamiltonian formulation of [2| presented an 

n n 

obstruction to the Jacobi identity [34|. This term was replaced in [12] in order remove this 
problem, but it was then reported in [3] that the new term also presents an obstruction, 



viz. that given by flA22j) below. One can rescue the Hamiltonian theory by requiring all 
functionals to depend on fields B such that V • B = 0, but because the orginal program 
begun in was to construct truly gauge-free field theories in terms of noncanonical 

Poisson brackets, this taint was a disappointment. The same obstruction appeared in the 
context of magnetohydrodynamics [2], but a way to remove it was obtained in [3]. To date, 



the best fix for the taint of the Vlasov-Maxwell bracket is given in 



251 ] by using Dirac 
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constraint theory, which replaces V • B = by a boundary condition at infinity. 

Appendix I A II contains the details of the onerous calculation first performed by the author 
in 1981 (reported in detail here for the first time, as it was originally done), leading to the 
result of (IA22|) that appeared in . This is followed in Appendices IA2I and IA 31 by a direct 
proof of the brackets for the spin and monopole Vlasov-Maxwell theories, respectively. 



1. Jacobi identity for the Vlasov-Maxwell bracket 

For convenience the charge, mass, and a factor of Arc are scaled out to obtain the Vlasov- 
Maxwell bracket for the fields f(z,t), -E(x, t), and -B(x, t) in the following form: 

{F, G} = j d 6 zf[F f , G f ) c + f[F f , G f ] B + f [G E ■ d v F f - F E ■ d v G f ) 

+ [d 3 xF E -V x Gb-Ge-V x F B 



=: {F, G} c + {F, G}b + {F, G} Ef + {F, G} EB , (Al) 

where F f := 6F/5f, [f,g) c : = V/ • d v g - Vg ■ d v f, d v := d/dv, [f,g] B := B ■ (d v f x d v g), 
and {F, G} c etc. are obvious from context. Also, note, boldface has been removed since the 
formulas are busy enough. Since charge can be scaled out in this manner, it is evident that 



the validity of the Jacobi identity is independent of the sign of t 



re species charge. 



Above, the term {F, G} B , the Marsden-Weinstein term [12J, has been separated out 
because it will be seen to be the source of the failure of Jacobi identity unless V • B = 0. 
Considering the combination {F, G} c + {F, G}b together would simplify the calculation 
somewhat. 

The Jacobi identity is 

{{F, G}, H} :=: {{F, G} c , H} + {{F, G} B , H} + {{F, G} Ef , H} + {{F, G} EB , H} 

= {{F, G} c , H} c + {{F, G} c , H} B + {{F, G} c , H} Ef + {{F, G} c , H} EB 

V v ' V v ' S v ' V V ' 

12 3 4 

+ {{F, G} B , H} c + {{F, G} B , H} b + {{F, G} B , H} Ef + {{F, G} B , H} EB 

S V ' V v ' * v ' V v ' 

5 6 7 8 

+ {{F, G} E f, H} c + {{F, G} Ef , H} B + {{F, G} Ef , H} Ef + {{F, G} Ef , H} EB 
* v ' > „ ' > „ ' > v ' 

9 10 11 12 

+ {{F, G} EB , H} c + {{F, G} EB , H} B + {{F, G} EB , H} Ef + {{F, G} EB , H} EB 

S v ' V v ' ^ v ' N v ' 

13 14 15 16 
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where the symbol :=: means modulo cyclic permutation; 'eye' will be used to denote cyclic 
permutation when only one side is to be permuted. 

To prove the Jacobi identity for a generic Poisson bracket of the form {F, G} = [F^J'G^) , 
with cosymplectic operator J ', one must calculate the functional derivative 5{F, G}/5ip- This 
derivative has three contributions: one from F^, one from J ', and one from G^. The first 
and last give rise to second functional derivatives. Proofs of the Jacobi identity are greatly 
simplified by the following: 



Bracket Theorem (|3() To prove the Jacobi identity for generic brackets of the form 
{F, G} = (F^IJ'Gtp) one need only consider the explicit dependence of J on ip when taking 
the functional derivative 5{F, G}/S^. 

Proof. The formal proof uses the anti-self-adjointness of J and the self-adjointness of 
the second functional derivative. With these symmetries it can be shown that all second 
functional derivative terms cancel. □ 

In what follows 5{F, G}/5ifj = ... denotes the functional derivative modulo the second 
derivative terms. Equation ( lAlj) gives 



S{F,G} e ^ 6{F,G} C ^ S{F,G} C ^ 

—Sj—-[Ff,G,]e, -0, — ^--0, (A2) 

6{F,G} B . rc , „ , 6{F,G} B . r 5{F,G} B 



Sf 



[Ff,G f ) B , 6 -^j^=Jdvf(d v F f xd v G f ), ^P = 0, (A3) 

6 A£^hl = {GE . dvFf _ FE . dvGfh ^§^ = 0, ^§^ = 0, (A4) 

6{F, G} EB ^ 5{F, G} EB ^ n 6{F, G} EB ^ 

If -°' SB -°' 6E "°- (A5) 

The following are immediate: 

• Term 1 vanishes because {F, G} c is Lie- Poisson; i.e., using the first of (1A2j) . 5{F, G} c /5f = [Ff, Gf], 
which, when inserted into {{F, G} c , H} c and cyclicly permuting, vanishes by virtue of 

the Jacobi identity of [ , ] c . 

• Term 4 vanishes by the Bracket Theorem because of the second two equations of (1A2j) . 

• Term 12 vanishes by the Bracket Theorem because of the second two equations of f ]A4j) 

• Terms 13-16 vanish by the Bracket Theorem because {F, G}eb has no explicit depen- 



dence on /, E, or B, i.e. because of flA5j) . 
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Remark. One can organize Jacobi identity calculations at the outset by grouping together 
all like terms that can possibly cancel. For example terms with the same functional deriva- 
tives of F, G, and H must be considered together. Sometimes other considerations can aid 
in the grouping of terms. When terms are grouped appropriately, failure of a class of terms 
to cancel is a proof of the failure of Jacobi. In the heading below /// means that only 
function derivatives with respect to / occur, etc. 

Term 6 (///) 
Using the first equation of (1A3j) gives 



6 := 



=: J f[B-(d v F f xd v G f ),H f } B 
J fB- (d v (B ■ (d v Fj x d v G f ) ) x d v H } 



f BiB r 6ijk€ rs t 



d fdFfdGA dH 



dvj \ dv s dv t 



fBiB r eijke 



ijk^rst 



dv k 
+ 



d 2 F f dGfdH f dF f d 2 G f dH f 



1 



f BiB r €ijk€ r 



St 



dvjdv s dv t dvk dv s dvjdv t dvk 
d 2 F f dG f dH f dHf d 2 F f dGj 



dvjdv s dv t dvk dvk dv s dvj dv t 



(A6) 
(A7) 



where ( 1A7I) follows from ( 1A6I) by permuting the second term of ( 1A6I) . shifting the indices 
according to s — Y k, k — >■ t, t — > j, j — >■ s, and i -f-)- r, and using the antisymmetry of the 
Levi-Civita symbol. 

Note, the above procedure is common in this game and of general utility, so it is recorded 
in the following: 

Lemma 1 If two terms can be made to cancel by permuting one of them, then all terms 
cancel. 

Proof. By writing out all six terms by permuting FGH — > GHF — > HFG, one observes 
they cancel in pairs. □ 

Remark. Term 6 vanishes without any assumptions on B, i.e. V ■ B = is not required. 

Term 11 (EEf) 
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Using the first and last equations of (IA4I) yields 
11:=: ff Er 9,WG} " 



<9t> 5/ 

:=: J fH E -d v (G E ■ d v F f - F E ■ d v G f ) 

:=: J f (H E ■ d v {G E ■ d v F f ) - G E - d v {H E ■ d v F f )^} :=: 



(A8) 



where the last equality follows because com{G E ■ d v ,H E ■ d v ) = 0, where com means com- 
mutator. 

Terms 2 and 5 (///) 
Remark. Terms 2 and 5 have been grouped together because both give rise to terms in the 
Jacobi identity involving Ff, Gf, and Hf. 



Using dX2D in 2 + 5 :=: {{F, G} c , H} B + {{F, # } c gives 



2 + 5 :=: / /S • (^[F/, G f ] e x ^F/) + / [S • x d v Gf), G f ] c 



fB t e 



ijk 



dHf 



dF) 

dvi 



G f 



+ 



/%& ~ — B 



d ( u dF f dGf\dH. 



Bi- 



d fdFfdG f \ OH 



\ox£ \ ovj ovk J ovi ove \oVj ovk J oxg 



/ Bi €ij k 
+ / / Cijk 



dH 



f 



OBi fdF f dGfdHf 



dF 



+ 



7» 



+ 



<9t> ,• (9f fc ' 



(A9) 
(A10) 



Upon defining F = d v Ff, Q = d v Gf, and H = d v Hf, and using the Leibniz rule for [ 
Line (IA9D can be rewritten as 



J fB^^i-HklF^Gfl + nklF^g^ + gklF^Hfl + FAGk^Hfl) :=:0, (All) 

where upon permutation the first and fourth terms cancel, as do the second and third. Thus, 

OBi fdFfdGfdHf 



2 + 5:=: / ft 



dxi \ dvj dvk dve 



(A12) 



From (|A12j) . the only remaining term is from Line (lAlOj) . This term can be rearranged to 
yield 

2 + 5 + eye = J fV-B [(d v Ff x d v G f ) ■ d v H f ) , (A13) 
which is a consequence of the following: 
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Lemma 2 For any three vectors J 7 , Q, and 7-L, and vector field B in M 3 , 

H -VB ■ (F x g) + T -VB ■ (G xH) + G -VB ■ {H x JF) = (V • B)(U ■ (J 7 x Q)) (A14) 
where U-VB -C = H&jdiBj. 

Remark. Terms 2 and 5 could have been combined with Term 6. They were considered 
separately to pinpoint, as will be seen, that they are the sole terms that violate the Jacobi 
identity without V • B = 0. 

Terms 7 and 10 (Eff) 
Inserting the first and last equations of (IA3j) and the first equation of (IA4j) into 7+10 :=: 
{{F, G} B , H} Ef + {{F, G} Ef , H} B gives 

7+10 :=: J fH E -d v (B-(d v F f xd v G f )) 

+fB- (d v {-G E ■ d v F f + F E ■ d v G f ) x d v H f j) 

'■='■ J f H Ei d v£ ( y B i e i jkd V jFfd vk G f) 

—fBi e,ijk (—d V j(G Ee d ve Ff)d vk Hf + d V j(F E£ d ve Gf)d vk Hf^ 

'■='■ J f Bi€ij k H El (d vk Gfd v£ d V jFf + d vj Ffd vi d vk Gf) (A15) 
—fBieijk(—G E ed V kHfd V jd ve Ff + F Ee d V kHfd V jd v eGf) . (A16) 

Now upon permutation, the first term of ( )A15j) is seen to cancel the second of ( ]A16j) and 
the second term of (1A15|) cancels the first term of (1A16j) . Thus 7 + 10 :=: 0. 



Terms 3, 8, and 9 (Eff) 
Using the first and last equations of (1A2j) and (lA4j) and the first and second equations of 
3]) in 3 + 9 gives 



3 + 9:=:-Jf(H E -d v [F f , G f ] c + [H E ■ d v G f , F f ] e - [H E ■ d v F f , G f ] c ) (A17) 
Using Lemma 3 below in flA17j) with C equal to H E gives 



3 + 9 :=: J f (V x H E ) ■ (d v F f x d v G f ) . (A18) 

Lemma 3 For any vector field C(x) and phase space functions f and g, 

C ■ d v [f, g) c = [C ■ d v f, g) c + [/, C ■ d v g] c + (V x C) ■ (d v g x d v f) . (A19) 
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Proof. With d vi := d/dvi and <9 ; 



xi 



d/dx. 



C-d v [f,g] c = C t ([d v J,g] c +[f,d vl g] c ) 

= [C ■ d v f, g] c - [d, g] c d v J +[f,C- d v g] c - [f, d] e d vi g 

= [C ■ d v f, g] c - {d vi f){d vj g)d xj Ci + [f,C ■ d v g} c + {d vi g){d vj f)d xj Ci 

= [C ■ d v f,g} c + [f,C- d v g] c + (d^C^dW^g) - (d vi g)(d vj f)) 

= [C ■ d v f, g] c + [f,C ■ d v g] c + (d vi g)(d vj f) (d xj Ci - d xi Cj) 

= [C ■ d v f, g] c + [f,C- d v g] c + (V x C) ■ (d v g x dj) . □ (A20) 



Equations flAT8l) and flA2T]) imply 3 + 8 + 9 + eye = 0. 

Remark. Observe the terms here, like Terms 7 and 10, are Eff terms. However, they have 
been grouped separately because there are 'other considerations' as mentioned above. The 
Terms 7 and 10 vanish with B, but the Terms 3, 8, and 9, do not. Thus terms of one kind 
cannot cancel terms of the other. 

Finally, from all of the above, the following is concluded: 

Main Theorem ({3]]) For the Vlasov-Maxwell bracket of $A1\) 



Remark. It is interesting to note that the other constraint, V • E = 4:irp } need not be 
satisfied for the Jacobi identity to hold. It turns out to be a Casimir invariant. 

2. Jacobi identity for the spin Vlasov-Maxwell bracket 

Writing {F, G} s vm = {F, G}vm + {F, G} s and using :=: as defined in Appendix IA II 



Now consider Term 8 



S 




(A21) 




(A22) 



{{F, G}svm, H} s vM {{F, G}vm, H}vm + \{F, G} s , H}vm + 



+{{F,G} VM ,H} S + {{F,G} S ,H} 



'■='■ {{F, G} s , H} VM + {{F, G} VM , H} 



s ) 



(A23) 
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where the second equality follows because of the Jacobi identity for Vlasov-Maxwell (assum 



ing solenoidal B) and the fact that {F, G} s is a Lie-Poisson bracket (see e.g. |4], |35j)- Thus 
it only remains to show that the cross terms cancel, which is facilitated by a the bracket 
theorem stated in Appendix viz., when functionally differentiating {F, G}vm and 
{F, G} s , which are needed when constructing the cross terms, one can ignore the second 
functional derivative terms. Using the symbol = again to denote equivalence modulo the 
second variation terms, 

5{F, G}vm 



Sf 
S{F,G} f 



= [F f , G f ] c + [F f , G f ] B + F E -d v G f -G E - d v F f (A24) 
= [F f ,G f ] s (A25) 



5f 

while all other needed functional derivatives vanish. Thus 



{{F,G} VM ,H} S :=: jd'z (/ [[F f ,G f ) c + [F f ,G f ) B ,H f ] s 

+f[F E ■ d v G f - G E ■ d v F f ,H f ] s ^ (A26) 
{{F,G} S ,H} VM :=: jd 9 z (f[[F f ,G f ] s ,H f ] c + f[[F f ,G f ] s ,H f ] 



+fH E -d v [F f ,G f ] s ) (A27) 
The first lines of ( 1A26j) and (1A27I) cancel by virtue of the Jacobi identities for the brackets 



■ yE ,s on functions, while the second line of ( ]A26j) cancels upon permutation of the second 



term. Similarly, the second term of (1A27I) vanishes. 

3. Jacobi identity for the monopole Vlasov-Maxwell bracket 

For this case the Gaussian units of the text are used, i.e. the factors of Air are reinserted 
and both the usual and monople charges are manifest. Let {F, G} m = {F, G}vm + {F, G}m, 
where 

{F,G} M = fd 6 zf[F f ,G f } E + ^f(G B -d v F f -F B -d v G f ), (A28) 
and [ , ]e is defined by fl79|) . Thus, the Jacobi identity has four terms to consider 
{{F, G} m , H} m :=: {{F, G} V m, H} m + {{F, G} M , H} m 

= {{F, G}vm, H} VM + {{F, G} M , H} vm 



1 2 



{{F, G} VM , H} M + {{F, G} M , H} M , (A29) 



1 
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where the symbol :=: is defined in Appendix IA II 

Term 1 

From the [3| (cf. Appendix IA II) 



{{F, G} VM , H} VM + C yc = ^ [d 6 zfV-B [(d v F f x d v G f ) ■ d v H f ] . 



(A30) 



As in Appendix IA II 5{F,G}/5ip = ... denotes the functional derivative modulo the 
second derivative terms. The following will be needed: 



5{F, G}vm = + Gf]B + 4vre ^ Q ^ _ ^ ^ 



m 



4irg 



= [F f , Gj]e H (G B ■ d v Ff - F B ■ d v Gf) 

m 



8{F,G} M 
Sf 

5{F, G}vm 

Fd = ~2 I a ' V J ° vP f X ° vU f ' c / 

5B m 2 J SE 

S{F,G} M , g p „^ 5{F,G} M 



d 3 vfd v F f x 0„G 



VM 



— ^ [d 3 vfd v F f xd v G f , 
m 2 J 





: 0. 



(A31) 
(A32) 
(A33) 

, , , ,.. . . , . „ . (A34) 

5E m 2 J oB 

The Poisson bracket that generates the l v x' part of the generalized Lorentz force is 
[f, 9] = [f, g]c + [f, 9\b + [f, 9]e- Because of 

[[f,g]E,h]E + cyc = and [[/, g} E , h] B + [[/, g] B , h] E + eye = . (A35) 
the following holds: 

[[/, g],h} + eye = [[/, (?] c , /i] B + [[/, ^] B , /i] c + [[/, g] c , h] E + [[/, g] E , h] c + eye 

= (eV ■ B — gV ■ E) [(d v f x 3 v g) ■ d v h] /m 2 (A36) 

The first term of the above is the source of the RHS of ( 1A30I) . 
Now consider the remaining terms of (1A29j) . 

Term 4 

Equations flA3~2]) and flA28l) give 



{{F,G} M ,H} 



d 6 zf 



[Ff, Gf} E H (G E ■ d v Ff — F B ■ d v Gf) , H 



m 



+ ^lfH B -d v ([F f ,G f } E + ^(G B -d v F f -F B -d v G f )) . (A37) 
m \ m J 

Equation AA37J) has three kinds of terms. Consider first the ///-terms: 

J d 6 zf [[F f ,G f ] E ,H f ] E = ^ J d 6 z f E ■ (d v (E ■ (d v Ff x d v Gf) ) x d v Hfj 

= J d e z f eijke rst EiE r d vk Hf (d V jd vs Ffd vt Gf 

+d vj d vt G f d Vs F f ) (A38) 
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Upon permutation and reindexing, the two terms of (1A38j) cancel. This, of course, follows 
immediately from flA35j) - the above serves as a proof that [[f,g]E,h]E + eye = 0. Next 
consider the BBf-term: 

16ttV 



rrr 



J d 6 z fH B ■ d v (G B ■ d v F f - F B ■ d v G f ) 



(A39) 



This vanishes upon permutation because H B • d v and F B ■ d v commute. Now all that remains 
of Term 4 is the Bf /-term: 

Aug 



m 



J d 6 zf [(G B ■ d v F s - F B ■ d v G f ,H f ] E + f H B ■ d v [F f ,G f ] E 



(A40) 



This term is of the same form as the EFF term of the MV-bracket (terms 7 and 10), and 
vanishes for the same reason. Therefore, Term 4 vanishes. 
Now consider the two cross terms. 



Term 2 is 

{{F,G} M ,H} VM :=: J d 6 zf 

+ f 



Terms 2 and 3 



4ng 



[Ffi Gf]E + — (G B ■ d v Ff — F B ■ d v Gf) , Hf 



(A41) 
(A42) 



[•^/j Gj]e + — (G B ■ d v Fj — F B ■ d v Gf) , Hf 
+ ^fH E -d v ^F f ,G f } E + ^ L (G E • d v Ff — F E • d v Gf) \ (A43) 

(A44) 



- -fd v H f -^- I d 3 vfd v F f x d v G } 



- An ftV x H 



B 



d 3 v f d v F f x d v G 



while Term 3 is 



{{F,G} VM ,H} M :=: J d 6 z f [[F f ,G f } c + [F f ,G f ] B 



Ane 



m 



(G E ■ d v Ff — F E ■ d v Gf) , Hj 



+ -^■fH B -a ([F f ,G f ] e +[F f ,Gf] B 

(G E -d v F f -F E -d v Gf)} 



Ane 



- —f d v Hf ■ d 3 vfd v F f xd v G f . 
m m 2 J 

Upon comparing Terms 2 and 3 some cancellations are immediate. 



(A45) 

(A46) 
(A47) 
(A48) 
(A49) 
(A50) 
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Using (|A36|) . 1st term of flAlT]) + 1st term of (lA46|) gives - Jr(V • E) (d v F f x d v G f ) ■ 
d v H f 



Using (fA35|) . 1st term of (TA42I) + 2nd term of (!A46l) = 
Lines f lAiij) + (|A50|) =0 



The terms of (1A49I) vanish because Hb ■ d v and Ge • d v commute. Likewise the last 
two terms of ( 1A43I) 



Applying the following: 

Lemma For any vector field C(x) and phase space functions f and g, 

C ■ d v [f, g] c = [C ■ d v f, g] c + [f,C- d v g] c + m~\V x C) ■ (d v g x d v f) 



(A51) 



which is not difficult to prove, to the last terms of (IA41I) + Line flA45j) + the last first term 
of (IA48D =0. 

There are six remaining terms. The last two terms of ( 1A42I) cancel the last term of ( 1A48I ). 
and the first term of ( 1A43I) cancels the two terms of (1A47I) . Thus in this case, the obstruction 
becomes 

{{F, G} mV M, H} mVM + eye = (A52) 

\ [d 6 zf {eV-B-gV- E) (d v F f x d v G f ) ■ d v H f , 

and it is concluded that the Jacobi identity still requires a solenoid constraint on eB — gE. 

Upon transforming to new variables e~E = eE+gB and eB = eB — gE, where e 2 = e 2 +g 2 , 
reproduces the Poisson bracket for Vlasov-Maxwell theory, which is possible for this single 
species case of Dirac's theory. In retrospect, the existence of this transformation precludes 
the necessity for the proof of the Jacobi identity; however, consistent with the goal of this 
entire appendix, viz. to demonstrate techniques of general utility rather than to present the 
most efficient proofs, we retain it here. 
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